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Although the link between long Gamma Ray Bursts (GRBs) and supernovae (SNe) 
has been established 1-4 , hitherto there have been no observations of the beginning of 
a supernova explosion and its intimate link to a GRB. In particular, we do not know 
however how a GRB jet emerges from the star surface nor how a GRB progenitor 
explodes. Here we report on observations of the close GRB060218 5 and its connection 
to SN2006aj 6 . In addition to the classical non-thermal emission, GRB060218 shows a 
thermal component in its X ray spectrum, which cools and shifts into the optical/UV 
band as time passes. We interpret these features as arising from the break out of 
a shock driven by a mildly relativistic shell into the dense wind surrounding the 
progenitor 7 . Our observations allow us for the first time to catch a SN in the act of 
exploding, to directly observe the shock break-out and to provide strong evidence that 
the GRB progenitor was a Wolf-Rayet star. 

GRB060218 was detected with the BAT instrument 8 onboard Swift 9 on 2006 February 18.149 
UT 5 . The burst profile is unusually long with a T 90 (the time interval containing 90% of the flux) of 
2100 ± 100 s (Fig. 1). The flux slowly rose to the peak at 431 ± 60 s (90% containment, times are 
measured from the BAT trigger time). Swift slewed autonomously to the newly discovered burst. 
The X-ray Telescope (XRT 12 ) found a bright source, which rose smoothly to a peak of ~ 100 counts 
s -1 (0.3-10 keV) at 985 ± 15 s. The X-ray flux then decayed exponentially with an e— folding 
time of 2100 ± 50 s, followed around 10 ks by a shallower power-law decay similar to that seen in 
typical GRB afterglows 14-15 (Fig. 2). The UltraVio let /Optical Telescope (UVOT 13 ) found emission 
steadily brightening by a factor of 5-10 after the first detection, peaking in a broad plateau first in 
the UV (31.3 ± 1.8 ks at 188 nm) and later in the optical (39.6 ± 2.5 ks at 439 nm). The light curves 
reached a minimum at about 200 ks, after which the UV light curves remained constant while a 
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rebrightening is seen in the optical bands, peaking again at about 700-800 ks (Fig. 2). 

Soon after the Swift discovery, low-resolution spectra of the optical afterglow and host galaxy 
revealed strong emission lines at a redshift of z — 0.033 (ref. 16). Spectroscopic indications of the 
presence of a rising supernova (designated SN2006aj) were found three days after the burst 6,17 with 
broad emission features consistent with a type Ic SN (due to a lack of hydrogen and helium lines). 

The Swift instruments provided valuable spectral information. The high energy spectra soften 
with time and can be fit with (cut-off) power-laws. This power-law component can be ascribed 
to the usual GRB jet and afterglow. The most striking feature, however, is the presence of a soft 
component in the X-ray spectrum, that is present in the XRT up to ~ 10, 000 s. The blackbody 
component shows a marginally decreasing temperature (kT ~ 0.17 keV, where k is the Boltzmann 
constant), and a clear increase in luminosity with time, corresponding to an increase in the apparent 
emission radius from R% B = (5.2±0.5) x 10 11 cm to R^ B = (1.2±0.1) x 10 12 cm (Fig. 3). During the 
rapid decay (t ~ 7000 s), a blackbody component is still present in the data with a marginally cooler 
temperature (kT = 0.10 ± 0.05 keV) and a comparable emission radius [R$ B = (6.5i^ 4 4 ) x 10 11 
cm]. In the optical/UV band at 9 hours (32 ks) the blackbody peak is still above the UVOT energy 
range. At 120 ks the peak of the blackbody emission is within the UVOT passband, and the inferred 
temperature and radius are kT = 3.7tg;g eV and R BB = 3.29to!93 x 10 14 cm, implying an expansion 
speed of (2.7 ± 0.8) x 10 9 cm s -1 . This estimate is consistent with what we would expect for a SN 
and it is also consistent with the line broadening observed in the optical spectra. 

The thermal components are the key to interpreting this anomalous GRB. The high 
temperature (two million degrees) of the thermal X-ray component suggests that the radiation 
is emitted by a shock-heated plasma. The characteristic radius of the emitting region, 
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-Rsheii ~ (i?/aT 4 ) 1//3 ~ 5 x 10 12 cm (E is the GRB isotropic energy and a is the radiation density 
constant), corresponds to the radius of a blue supergiant progenitor. However, the lack of hydrogen 
lines in the SN spectrum suggests a much more compact source. The large emission radius may be 
explained in this case by the existence of a massive stellar wind surrounding the progenitor, as is 
common for Wolf-Rayet stars. The thermal radiation is observed once the shock driven into the 
wind reaches a radius, ~ -R s heii, where the wind becomes optically thin. 

The characteristic variability time is -R s heii/ C ~ 200 s, consistent with the smoothness of the 
X-ray pulse and the rapid thermal X-ray flux decrease at the end of the pulse. We interpret this 
as providing, for the first time, a direct measurement of the shock break-out 19,20 of the stellar 
envelope and the stellar wind (first investigated by Colgate 21 ). The fact that -R s hcii is larger than 
Rbb suggests that the shock expands in a non-spherical manner, reaching different points on the 
-Rsheii sphere at different times. This may be due to a non-spherical explosion (e.g. the presence of a 
jet), or a non-spherical wind 22 ' 23 . In addition, the shock break-out interpretation provides us with a 
delay between the SN explosion and the GRB start to be £4 ks 24 (Fig. 1). 

As the shock propagates into the wind, it compresses the wind plasma into a thin shell. 
The mass of this shell may be inferred from the requirement that its optical depth be close 
to unity, M shcll 4 7r i? 2 hcll //t 5 x 1O~ 7 M (k « 0.34 cm 2 g -1 is the opacity). This implies 
that the wind mass-loss rate is M m M s h e ii t> w ind/-Rsheii ~ 3 x lO~ 4 M yr _1 , for a wind velocity 
"^wind = 10 s cms -1 , typical for Wolf-Rayet stars. Since the thermal energy density behind a 
radiation-dominated shock is aT 4 ~ 3 pv 2 (p is the wind density at R S \ 1C \\ and v s the shock velocity) 
we have p ~ 10 -12 gcm -3 , which implies that the shock must be (mildly) relativistic, v s ~ c. This 
is similar to GRB980425/SN1998bw, where the ejection of a mildly relativistic shell with energy of 
~ 5 x 10 49 erg is believed to have powered radio 25-27 and X-ray emission 7 . 
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The optical-UV emission observed at early time, t £ 10 4 s, may be accounted for as the 
low-energy tail of the thermal X-ray emission produced by the (radiation) shock driven into the 
wind. At later time, the optical-UV emission is well above that expected from the (collisionless) 
shock driven into the wind. This emission is most likely due to the expanding envelope of the star, 
which was heated by the shock passage to a much higher temperature. Initially, this envelope is 
hidden by the wind. As the star and wind expand, the photosphere propagates inward, revealing 
shocked stellar plasma. As the star expands, the radiation temperature decreases and the apparent 
radius increases (Fig. 3). The radius inferred at the peak of the UV emission, ~3x 10 14 cm, 
implies that emission is arising from the outer ~ An (R^) 2 / k ~ 10~ 3 M Q shell at the edge of the 
shocked star. As the photosphere rapidly cools, this component of the emission fades. The UV 
light continues to plummet as cooler temperatures allow elements to recombine and line blanketing 
to set in, while radioactive decay causes the optical light to begin rising to the primary maximum 
normally seen in SN light curves (Fig. 2). 

Since the wind shell is clearly larger than the progenitor radius, we infer that the star radius 
is definitely smaller than 5 x 10 12 cm. Assuming a linear expansion at the beginning (due to light 
travel time effects) one can estimate a star radius of i? s tar ~ (4 ± 1) x 10 11 cm. This is smaller than 
the radius of the progenitors of type II SNe, like blue supergiants (4 x 10 12 cm for SN1987A, ref. 
28) or red supergiants (3 x 10 13 cm). Our results unambiguously indicate that the progenitor of 
GRB060218/SN2006aj was a compact massive star, most likely a Wolf-Rayet star. 
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Fig. 1.— Early Swift light curve of GRB060218. GRB060218 was discovered by the BAT when 
it came into the BAT field of view during a pre-planned slew. There is no emission at the GRB 
location up to -3509 s. Swift slewed again to the burst position and the XRT and UVOT began 
observing GRB060218 159 s later. The BAT light curve is shown with green open squares. The XRT 
data are shown with black open circles. For each BAT point we converted the observed count rate to 
flux (15-150 keV band) using the observed spectra. Combined BAT and XRT spectra were fit with 
cut-off power-law plus a blackbody, absorbed by interstellar matter in our Galaxy [column density of 
(0.9 — 1.1) x 10 21 cm -2 ] and in the host galaxy at redshift z = 0.033. The host galaxy column density 
is Nfj = 5.0 x 10 21 cm -2 . Errors are at 1 a significance. At a redshift z = 0.033 (corresponding to 
a distance of 145 Mpc with H = 70 km s _1 Mpc -1 ) the isotropic equivalent energy, extrapolated 
to the 1-10,000 keV rest frame energy band, is Ei SO = (6.2 ± 0.3) x 10 49 erg. The peak energy in 
the GRB spectrum is at E p = 4.9l ; 3 keV. These values are consistent with the Amati correlation, 
suggesting that GRB060218 is not an off-axis event 10 . This conclusion is also supported by the lack 
of achromatic rise behavior of the light curve in the three Swift observation bands. The BAT fluence 
is dominated by soft X-ray photons and this burst can be classified as an X-ray flash 11 . 
A V band light curve is shown with red filled circles. For clarity the V flux has been multiplied 
by a factor of 100. Magnitudes have been converted to fluxes using standard UVOT zeropoints and 
multiplying the specific flux by the filter Full Width at Half Maximum. Gaps in the light curve are 
due to the automated periodic change of filters during the first observation of the GRB. 
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Fig. 2. — Long-term Swift light curve of GRB060218. Upper panel: the XRT light curve 
(0.3-10 keV) is shown with open black circles. Count rate-to-flux conversion factors were derived 
from time-dependent spectral analysis. We also plot with open black squares the contribution to the 
0.3-10 keV flux by the blackbody component. Its percentage contribution is increasing with time, 
becoming dominant at the end of the exponential decay. The X-ray light curve has a long, slow 
power-law rise followed by an exponential (or steep power-law) decay. At about 10,000 s the light 
curve breaks to a shallower power-law decay with index —1.2 ±0.1 characteristic of typical GRB 
afterglows. This classical afterglow can be naturally accounted for by a shock driven into the wind 
by a shell with kinetic energy E s ^ e \\ ~ 10 49 erg. The t _1 flux decline is valid at the stage where 
the shell is being decelerated by the wind with the deceleration phase beginning at tdec ~ 10 4 s for 
M ^ 10 _4 (w w i n d/10 8 cms -1 ) Moyr^ 1 , consistent with the mass-loss rate inferred from the thermal 
X-ray component. 

Lower panel: the UVOT light curve. Filled circles of different colors represent different UVOT filters: 
red - V (centered at 544 nm); green - B (439 nm), blue - U (345 nm), light blue - UVW1 (251 
nm); magenta - UVM1 (217 nm) and yellow - UVW2 (188 nm). Specific fluxes have been multiplied 
by their FWHM widths (75, 98, 88, 70, 51 and 76 nm, respectively). Data have been rebinned to 
increase the signal to noise ratio. The UV band light curve peaks at about 30 ks due to the shock 
break-out from the outer stellar surface and the surrounding dense stellar wind, while the optical 
band peaks at about 800 ks due to radioactive heating in the SN ejecta. 
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Fig. 3. — Evolution of the soft thermal component temperature and radius. Upper panel: 
evolution of the temperature of the soft thermal component. The joint BAT and XRT spectrum 
has been fit with a blackbody component plus a (cut-off) power-law in the first ~ 3, 000 s (see also 
the caption of Fig. 1). The last point (circled in green) comes from a fit to the six UVOT filters, 
assuming a blackbody model with Galactic reddening [E(B — V) = 0.14] and host galaxy reddening. 
This reddening has been determined by fitting the Rayleigh- Jeans tail of the blackbody emission at 
32 ks (9 hours). The data require an intrinsic E(B — V) = 0.20 ± 0.03 (assuming a Small Magellanic 
Cloud reddening law 18 ). Lower panel: evolution of the radius of the soft thermal component. The 
last point (circled in green) comes from the fitting of UVOT data. The continuous line represents a 
linear fit to the data. 
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